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Abstract 

Parton distributions in the small x region are numerically predicted by using 
a modified DGLAP equation with the GRV-like input distributions. We find that 
gluon recombination at twist-4 level obviously suppresses the rapid growth of parton 
densities with x decrease. We show that before the saturation scale is reached, 
saturation and partial saturation appear in the small x behavior of parton distribu- 
tions in nucleus and free proton, respectively. The antishadowing contributions to 
the saturation phenomena are also discussed. 



PACS numbers: 12.38.Bx, 24.85. +p 



1 Introduction 

The QCD evolution equations for parton densities at twist-2 level, both the DGLAP 
equation [1] and BFKL equation [2] predict a rapid increase of the parton densities in 
the small x region due to parton splitting, and the unitarity limit is violated. Therefore, 
the corrections of the higher order QCD effects, which suppress or shadow the growth 
of parton densities, become a focus of intensive study in recent years. An important 
character at the small x limit is that the suppressed parton distributions approach a 
limiting form and unitarity is restored. This is called saturation. 

There are various ways to define and analyze the saturation phenomena based on per- 
turbative QCD [3,4]. The shadowing corrections of gluon recombination to the integrated 
parton distributions were mainly studied by adding nonlinear terms in the DGLAP evo- 
lution equation in the coUinear factorization scheme. A pioneering work in this aspect 
was derived by Gribov, Levin and Ryskin in [5] and by Mueller and Qiu [6] at the twist-4 
level. The GLR-MQ equation sums the contributions of gluon recombination diagrams 
using the AGK (Abramovsky, Gribov, Kancheh) cutting rule [7]. In the next step, the 
contributions of multi-parton correlations are summed by using the Glauber model in 
Mueller's works [8] and this Glauber-Mueller equation reduces to the GLR-MQ equation 
at the twist-4 approximation. 

Recently, the predictions of the GLR-MQ equation for the gluon saturation scale were 
studied in [9] . However, the apphcations of the AGK cutting rule in the GLR-MQ equation 
was argued in a more general consideration by two of us (WZ and JHR) in [10], where the 
Feynman diagrams are summed in a quantum field theory framework instead of the AGK 
cutting rule. We shall refer to this evolution equation as the MD-DGLAP equation. A 
major difference among the above mentioned nonlinear equations is that the momentum 
conservation is restored in the MD-DGLAP equation by the antishadowing corrections. 



which may change the predictions of the GLR-MQ equations. 

The purpose of this work is to study the behavior of the parton (quark and gluon) 
distributions in and x at high gluon density using the MD-DGLAP equation. It is 
know that the solutions of QCD evolution equations are sensitive to the input parton 
distributions. Works [9] use the CTEQ input distributions [11] at Q\ = lAGeV'^ to fit 
the HERA data for the structure function ^2(2;, Q^) of proton. Then they predict the 
saturation scales in proton and nucleus by evolving backward from a higher Q'^ scale, 
where the nonlinear terms in the equation can be neglected. Therefore, one can linearly 
add the input distributions of the nucleons in a nuclear target. However, the backward 
evolution paths for gluon and sea quarks are not unique and hence there are uncertainties 
in the results. Different from works [9], we use the GRV model [12] for input distributions, 
where the evolution begins in a very low Qq<1 GeV^. All partons in the GRV input 
take the valence-like form and it implies the finiteness of parton number and low density 
of partons. Therefore, we can construct the input distributions in the nucleus using the 
input set for proton and evolve them according to the standard evolution technique. 

We fit the parameters in the input distributions using the HERA data [13] for both 
^2(3;, (5^) and dF2{x,Q^)/dlii.Q'^ in proton in a limiting kinematical region. Then we 
predict the small x behaviors of parton distributions at different scale in proton and 
nucleus, respectively. We find obvious screening effects in quark and gluon distributions. 
We also show that a partial saturation appears in the parton distributions of proton. In 
particular, a flatter plateau appears in the region of smaller x and lower in middle and 
heavy nuclei. However, we have not found the saturation scale in the expected domain 
according to the definition in literatures [5,14]. We introduce the new scales Qj^ and Xg 
to describe the small x behaviors of parton distributions in the leading recombination 
region. We also study the contributions of the antishadowing terms in the MD-DGLAP 



equation. Antishadowing compensates the lost momentum in shadowing. Although this 
lost momentum is little, our calculations show that the contributions of the antishadowing 
terms to the saturation phenomena can not be neglected. 

The paper is organized as follows. In section 2 we introduce the MD-DGLAP equation 
and compare it with the GLR-MQ equation. In section 3 we fit the parameters in the 
GRV-like input distributions in the proton using the MD-DGLAP equation and HERA 
data, then we predict the parton distributions beyond the HERA region. In section 4 we 
discuss the parton distributions in the nucleus. The discussions and a summary are given 
in section 5. 



2 The MD-DGLAP equation 

As we know that the DGLAP equation produces a rapid growth of gluon density toward 
smaller values of x. The gluons therefore must begin to spatially overlap and recombine 
in a thin target disc. In works [10], the corrections of parton recombination to the QCD 
evolution equation are considered by summing up all possible twist-4 cut diagrams in the 
LLA((5^). In the derivation of the equation, the time-ordered perturbation theory instead 
of the AGK cutting rule is used to pick up the contributions of the leading recombination 
diagrams. In consequence, the corrections of the gluon recombination to the evolution of 
parton distributions with Q"^ are described by the following modified DGLAP equation 
[10] 
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for sea quark distributions, where P^^ are the evolution kernels of the linear DGLAP 
equation and the recombination functions 
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The nonlinear coefficient K in Eq. (1) depends on the definition of double parton 
distribution and the geometric distributions of partons inside target. For simplicity, we 
take as a free parameter in this work. Comparing with the GLR-MQ equation [6]: 
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there are following properties in Eq. (1): 

(i) The third term on the right-hand side of Eq. (1) is positive and it is called as the 
antishadowing, while the negative fourth term arises from the shadowing correction. The 
coexistence of shadowing and antishadowing in the QCD evolution of the parton densities 
is a general requisition of the local momentum conservation. We emphasize that the 



shadowing and antishadowing terms are defined on different kinematics domains [x, 1/2] 
and [x/2,x], respectively. Thus, the net recombination effects in Eq. (1) are not only 
related to the value of gluon density, but also depend on the slope of the gluon distribution 
in the space [x/2, x], i.e., a steeper (or flatter) gluon distribution has a stronger (or weaker) 
antishadowing effect. On the other hand, the AGK cutting rule is used in the derivation 
of Eq. (2) and where it is assumed that both the positive and negative contributions of 
the Feynman diagrams have the same kinematics domain. This results in the violation of 
the momentum conservation in the GLR-MQ equation. 

(ii) The GLR-MQ equation (2) takes the double leading logarithmic approximation 
(DLLA) for both and 1/z, where one keeps only the ln((5^//x^) ln(l/2;) factor in the 
solutions of the evolution equation or, equivalently, takes only the terms having 1/z — 
xi/x factor. On the other hand, the MD-DGLAP equation (1) is derived under the 
leading logarithmic approximation for Q'^. Therefore, Eqs. (Ic) and (Id) contain the 
terms beyond the leading 1/z approximation. One can find that the contributions from 
these terms do not vanish even in the small x region, since z runs from a; to 1 as Xi varies, 
that is z is not restricted in a smaller-z region. 

(iii) Note that the sea quark evolution in Eqs. (1) and (2) take different forms. The 
reason is that the transition of gluon— > quarks is suppressed in the DLLA-manner. The 
DLLA diagram only contains the gluon ladders and any transitions of gluon quark break 
the gluon ladder-structure. Therefore, a special technique is used to include the corrections 
of gluon recombination to the quark distributions in the GLR-MQ equation. However, 
this extra technique is unnecessary in the derivation of the MD-DGLAP equation, since 
we can produce the evolution equations for gluon and sea quarks in a same framework at 
the LLA{Q^). 

Therefore, we use the MD-DGLAP equation to study the saturation phenomena in- 



stead of the GLR-MQ dynamics. We shall show the different predictions of the evolution 
equations with or without antishadowing terms in Sec. 5. 



3 The solutions of the MD-DGLAP equation in 

proton 

The numerical solutions of evolution equation depend sensitively on the input parton 
(gluon and quark) distributions at a lower scale Qq. In principle, they are not calculable 
within perturbative QCD but are determined by data. Because the electromagnetic probe 
can not directly measure the gluon density, the input gluon distribution has a larger 
uncertainty. The early data for the DIS structure functions can always be fitted by using 
the linear DGLAP equation provided that a satisfying input gluon density is assumed. 
However, new HERA data in small x region, in particular the slopes of the structure 
function dF2/dlnQ'^ have restricted the above mentioned uncertainty. In fact, the global 
analysis of the HERA data using the DGLAP dynamics are given by the MRST [15] and 
CTEQ [11] collaborations. Where a MRST2001 set shows the negative value of gluon 
distribution in Q'^ < IGeV^. It means that if a positive input gluon distribution below 1 
GeV^ is taken, the screening corrections to the evolution equation are useful. 

In this work, we use the GRV-like input distributions. In the GRV model [12], the 
parton distributions are evolved from a very lower resolution scale (but larger than the 
QCD parameter Aqcd)- A specific assumption of the GRV model is that the input parton 
distributions take the simple valence-like distribution form. The GRV model with the 
linear DGLAP evolution equation gives a good description for proton structure function 
F2{x,Q'^) in a broad region but the fit in the slope of F2{x,Q'^) at lower remains to 
be improved (see the dashed curves in Fig. 2). We shall show that the corrections of the 
gluon recombination improved the fit. In practice, we use the initial valence quark- and 
gluon-densities in the GRV98LO set [16] as the input distributions at Ql = O.SiGeV^, 
i.e.. 
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and the mass of charm quark is 1.4 GeV. In the meantime, we let the parameters in the 
sea quarks distributions to be determined by the HERA data [13] in the MD-DGLAP 
evolution equation. The results are 



In this fit, the nonlinear coefficient K in the evolution equation is taken sls K — 0.0014. It 
implies that the nonlinear recombination corrections can not be neglected in the HERA 
data [13]. The results are shown in Figs. 1 and 2 (sohd curves), where the dashed lines 
are the fitting results of the linear DGLAP equation using the GRV98LO [16] as the input 
distributions. One can find that the contributions of gluon recombination improve the fit 
in Q'^ > IGeV"^. There are derivations between the fit and data in < IGeV"^, they 
imply that the corrections of beyond LLA((5^) are necessary at very lower Q^ . However, 
as we shall show that our main conclusions in this work are insensitive to the choices of 
parameters in the GRV-like input. 

Now we give the predictions of the gluon and quark distributions beyond HERA region 
in Figs. 3-6 (solid curves). The parton distributions from the DGLAP equation with 
GRV98LO are also plotted (dashed curves) for comparison. The results show obvious 



x{d{x, Ql) + u{x, Ql)) = 0.9x°-°^(l - a;)*^ "(l - 3.6^^ + 7.8a;), 



x{d{x, Ql) - u{x, Ql)) = 0.23x°-^^(l - x)^^-%l - 12.0^/x + 50.9x). 



(4) 



suppression in both the quark and gluon distributions and a flat tendency at the small 
X limit. Figures 5 and 6 give the Q^-dependence of parton distributions at fixed values 
of X. The dashed curves (the predictions of the DGLAP dynamics) in Figs. 3 and 4 



have the exponential form (~ expfyln 1/x]) and it violates the unitarity. However, the 
solid curves (the results of the MD-DGLAP equation) show that the growth of xG{x, Q^) 
and xS{x,Q'^) is slower than ln(l/a;) in 1 < < lOGeV'^ and x < 10~^. This partial 
saturation behavior is softer than the predicted result by the DGLAP equation and the 
BFKL equation (~ x~^, A > 0). 

In hadron-hadron cross section, such as p-p, the Froissart boundary [17] requires 



at the high s limit, where s is the center-of-mass energy squared and x — Q^/s. The 
high energy limit implies the small x approximation when is fixed to be a few GeV^. 
Therefore, although the gluon- and quark-distributions in proton do not saturate at small 
values of x, but their partial saturation-behavior satisfies the Froissart boundary in the 
perturbative QCD means. 

As we know that the saturation scale Ql{x), which indicates the saturation limit, is 
usually defined [5, 14] to be 
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or equivalently 



nonlinear terms of Eq.(l) 
W,r- J ^ \ J 



= 1, (7) 



linear terms of Eq.(l) 
which requires that the nonhnear recombination effect in the MD-DGLAP equation fully 
balances the linear splitting effect. Thus the saturation limit is reached. However, we 
have not found such saturation solution in Eq. (1). It means that the net shadowing 
(i.e., shadowing- antishadowing) effect in the leading recombination approximation is not 
large enough to cancel the increase of parton densities with increasing Q^. In other 
words, the higher-order recombination contributions would become significant and should 
be included in the evolution equation near the saturation limit. However, in this paper 
we focus on the range where the gluon recombination begins to work. Therefore, we 
introduce the recombination scale Q%{x) instead of Ql{x) as follows: 



Wi 



R 



nonlinear terms of Eq.{l) 



Ql = (^s{Ql{x)). (8) 



linear terms of Eq.{l) 
The Q\{x) for the gluon distribution in proton is shown in Fig. 7 (solid curve). It is 
interesting that the line has a corner near Qji{x) ~ 2GeV^. The anomalous behavior of 
the line at Qr{x) < 2GeV^ is the results of the antishadowing corrections. In fact, the 
relative stronger antishadowing effect at lower locally raises the gluon distribution. 
It shifts the flatter plateau toward a larger x- value. We can also see this effect in the 
behaviors of Xs in Figs. 3, 8 and 9 (see following sections). 

We can find that the gluon recombination obviously suppresses the evolution of parton 
densities in x < 10~^ in proton. According to the definition, one can see Q%{x) < Ql{x). 
The evolution of the parton distributions from Qjiix) to Ql{x) is a complicated process. 
We shall show that the contributions of the leading recombination are also important in 
the nuclear target near Qj^. 



4 The solutions of the MD-DGLAP equation in 

nucleus 

Nuclear target is an ideal laboratory for the research of the saturation phenomena, 
since the gluon recombination corrections are enhanced due to the correlation of gluons 
belonging to different nucleons at a same impact in the nuclear target. It is well known 
that the nuclear shadowing is a complicated phenomenon which has two different sources 
[18]: (i) it perturbatively originates from gluon recombination in the infinite momentum 
frame of nuclear target, or from multiple scattering in the target- rest frame; (ii) the non- 
perturbative nuclear effects. The former is expressed as the nonlinear QCD evolution 
equations, while the later relates to the structure of the input parton distributions. We 
have interest in the separate relations of nuclear saturation phenomena with the above 
mentioned two sources. Therefore, at first step, we neglect the non-perturbative nuclear 
effects. The gluon distribution (3) implies a small total number and very low density 
of gluons, where the recombination corrections are negligible. This conclusion can be 
confirmed from the following Fig. 8, which indicates that the gluon distributions with 
and without recombination corrections are similar near the evolution start point. Thus, 
we can predict the nuclear parton densities using the input distributions (3-4) and the 
MD-DGLAP equation, where the nonlinear terms are multiplied by A^^^. 

Figures 8-11 and 12-15 are the results corresponding to Figs. 3-6 but for Csi{A — 40) 
and Ph{A — 208), respectively. Our results show the plateau in the parton distributions 
inx < 10~^ and 1 < < lOGeV^. We note that in the DGLAP dynamics the generated 
gluon and sea quark distributions in small x region will speed up transfer from the valence 
form into the power form also through a plateau. However, this plateau exists in a very 
narrow Q^-window and it is unstable. The plateau in Figs. 8-9 and 12-13 imphes a 
saturation behavior in the parton distributions at the small x limit. We use Xg to indicate 



this saturation effect and is defined as 



dxGix^ =0 (9a) 

X 
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X 

The value of Xg is in windows 10"'' < Xs < 10"** and 1 < < lOGeV^ for the middle and 
heavy nuclear target. The results show that the antishadowing effect impedes the shift 
of Xg toward smaller value with increasing Q'^ in the range of < 2GeV'^. We can also 
find that the altitude of the plateau almost grows linearly with increasing in Fig. 16. 
This consists with the geometric scahng [3,4] . 

The values of for the nuclear gluon distribution are plotted in Fig. 7 (dashed and 
dotted curves) using Eq. (8). The behavior of curves near Qji{x) ~ 2GeV'^ is due to the 
antishadowing corrections, which locally raise the gluon distribution. 

Since the MD-DGLAP equation is based on the coUinear factorization scheme and 
its predictions for the parton distributions are universal and independent of the concrete 
process, the above mentioned saturation phenomena can be checked up in ultra-relativistic 
heavy ion collisions such as the rapidity distribution and centrality dependence of particle 
production. 

However, we still have not found the saturation solution Ql of Eq. (7) in the available 
nuclear target. Obviously, this conclusion is different from the work [9], which uses the 
GLR-MQ equation to backward evolve the parton distributions from the input distribu- 
tions of a free proton at very higher Q'^, where the nonlinear terms are negligible. They 
give g2 ^ 3 - 20GeV^ at 10"^ < x < 10"^ for Pb(^ = 208) [9]. One can understand 
the above mentioned difference as follows. Comparing the sohd curves of Figs. 3 and 12, 



we find the gluon density in the nucleus at higher is much smaller than that in the 
proton even if the recombination terms in the evolution equation do not play a role at such 
higher . The reason is that the parton distributions at higher Q'^ always remember the 
recombination effects in their evolution process. Therefore, a larger, but not true, input 
distribution in the nuclear target predicted a stronger recombination effect in work [9], 
which may balance the parton splitting effect in the evolution, and give a solution in Eq. 



Now we consider the corrections of the non-perturbative nuclear effects in the nuclear 
input parton distributions. The fact that the structure functions of bound and free nu- 
cleons are not equal has been discovered long ago and it is called the EMC effect [19]. 
However, its dynamics is still an open problem since the several mechanisms dominate 
the EMC effect. Although there are some models to describe the non-perturbative nu- 
clear shadowing relating to the GRV- input distributions [20,21], they still have a larger 
uncertainty, particularly in the nuclear gluon distribution. In this case, we take following 
simplified factors Rs and Rq to describe the contributions of the non-perturbative nuclear 
shadowing hke Ref. [22], 



(7). 
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for sea quarks and 
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for gluon, where 
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in which Xn — 0.05, = 0.017^"^/^, and Kg — 0.09 are parameterized using the data 
about the EMC effect [19]; while we take Kq = Ks for the moment. In the meantime, we 
neglect the EMC effects in x > 0.05 since we focus the behavior of parton distributions 
at very small x limit. Using Eqs. (10)-(13) as the input distributions we evolve the MD- 
DGLAP equation. Our results are shown in Figs. 17-19 for Ca(^ = 40) and Pb(^ = 208). 
Comparing these results with Eqs. 7-9, we find that the differences between them are 
small. 

Since we have not enough data to fix the value of Kg, we take the parton distributions 
for I'h{A = 208) as an example and change Kq from 0.18 to to test the sensitivity to 
such choices. We find that the results are insensitive to the above mentioned change since 
the gluon distribution Eq. (3) m. x < 0.05 is very small. 

Now we can conclude that the saturation- or partial saturation phenomena in the 
small X behavior of parton distributions are dominated by QCD dynamics rather than 
non-perturbative nuclear shadowing corrections in the input distributions. 



5 Discussions and summary 

Let us discuss the antishadowing contributions. Part of momentum is lost due to 
the negative shadowing terms in the GLR-MQ equation. Most works do not make an 
attempt to correct the evolution equation from the momentum conservation since the lost 
momentum is only a few percent of the total momentum. However, as one of us (WZ) 
has pointed out that the antishadowing contributions, which balance the lost momentum 
in the shadowing effect can not be neglected [23]. For illustrating this point, we set 
the antishadowing terms to zero in the MD-DGLAP equation (1) and use it to evolve 
the parton distributions. In this case, we re-fit the input distributions (4) and take 
K = 0.000285 using the HERA data [13]. New parton distributions in proton, Ca(A = 40) 
and Pb(A = 208) are plotted by dotted curves in Figs. 3-4, 8-9 and 12-13, respectively. 
There are following two distinguishing features as compared with solid curves: 

(1) The results show that the parton distributions with antishadowing effect are lower 
than that without antishadowing effect. One can understand the above-mentioned results 
as follows. The main HERA data, which are used to fit the nonhnear parameter K have 
shown the steeper form at > IGeV^ and x > 10~^, where the antishadowing correction 
is larger. Therefore, to fit the same HERA data, the value of K in the evolution equation 
with the antishadowing terms should be larger than that in the evolution equation without 
the antishadowing terms, since in the former case an obvious antishadowing effect partly 
cancels the shadowing effect. On the other hand, the net screening effect fiats the gluon 
distribution when x enters into a smaller x region and therefore weakens the antishadowing 
effect. Therefore, a more stronger net screening effect appeared in Eq. (1). Comparing 
the solid curves with the dotted curves in Figs. 3 and 4, one can understand that the 
antishadowing corrections make the growth of the solid curves slower than that of the 
dotted curves. 



The above mentioned differences between the effects with antishadowing and without 
antishadowing also appear in the predictions of the MD-DGLAP equation (with antishad- 
owing) and the GLR-MQ equation (without antishadowing). Using the same program as 
we used in the MD-DGLAP equation, we fit the relating parameters in the GLR-MQ 
equation when using the GRV-like input parton distributions. Then we predict the par- 
ton distributions in proton and Ph{A — 208) (sec Figs. 20 and 21). The results show 
that the screening effect in the GLR-MQ equation is weaker than the net-screening effect 
in the MD-DGLAP equation if we use similar input parton distributions. Although the 
plateau has appeared in the parton distributions of Ph{A — 208), we still have not found 
saturation solution in Eq. (6) for the GLR-MQ equation. 

(2) An other interesting difference between the predictions of the equations without 
and with antishadowing corrections is the position of Xg- In the former always 
moves toward smaller value with increasing (note that Xs runs out the diagrams, for 
example, at > 2GeV'^ in Fig. 12), while in the later case this shift is impeded by the 
antishadowing effect in the range of lower (Q^ < 2GeV'^). This phenomenon relates 
to Fig. 7, where the line has a corner near Qji{x) ~ 2GeV'^ due to the relative stronger 
antishadowing effect at lower locally raises the gluon distribution. 

As we know that the nuclear shadowing has different manifestations: the suppression 
in the usual or unintegrated parton distributions and their evolutions, the reduction of 
the structure functions and the screening effect in the cross sections. There are many 
different kinds of model to study the above mentioned shadowing phenomena. In this 
work we do not try to compare the modified DGLAP equations with other versions of the 
saturation model such as the JIMWLK equation [4] . The reason is that they have different 
research subjects: the former case discusses the shadowing in the parton distributions, 
which are independent of the process, while the later case treats the unintegrated parton 



distributions and the cross sections, which are process-dependent. 

Finally, we take the LLA((5^) in this work. The contributions beyond the LLA((5^) are 
necessary for the improvement of predictions in the range of very lower [Q"^ < IGeV^). 
However, it relates to the further study including the evolution dynamics at higher order. 

In summary, the parton distributions in small x region are numerically predicted by 
using a modified DGLAP equation with the GRV-like input distributions. We show that 
the gluon recombination at twist-4 level obviously suppresses the rapid growth of par- 
ton densities with x decrease. The growth of predicted parton distributions in proton 
towards small x is slower than ln(l/2;) in x < 10~^. In particular, a plateau is formed 
in the parton distributions of the nuclear target at small x limit and 1 < < lOGeV^. 
The altitude of the plateau almost grows hnearly with Q"^ increasing, i.e., the geometrical 
scaling. Thus, the parton distributions in proton and nucleus unitarize, and the Froissart 
boundary is not violated in the asymptotic regime of high density QCD. The saturation 
(or partial saturation) phenomena appear before the saturation scale Q^, where the gluon 
recombination correction fully balance the parton splitting effects. The predicted satu- 
ration phenomena can be checked up in ultra-relativistic heavy ion collisions such as the 
rapidity distribution and centrality dependence of particle production. 
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Figure Captions 

Fig. 1 The fits of the computed F2{x, Q^) in proton by the MD-DGLAP equation 
(solid curves) to (a) HI- and (b) ZEUS-data. The dashed curves are the DGLAP equation 
results from GRV98LO. 

Fig. 2 The fits of the computed dF2{x,Q'^)/dlnQ^ in proton by the MD-DGLAP 
equation (solid curves) to the HI data,. The dashed curves are the DGLAP equation 
results from GRV98LO. 

Fig. 3 The predictions for the gluon distribution function in proton. The solid (dotted) 
curves are the results of the MD-DGLAP equation with (without) antishadowing correc- 
tions; The dashed curves are the results of the DGLAP equation with the GRV98LO. 
Notice that the sohd and dashed curves have opposite concavities. 

Fig. 4 As Fig. 3 but for the sea-quark distribution function in proton. 

Fig. 5 The Q^-dependence of the gluon distribution function in proton at fixed values 
of X, the solid (dotted) curves correspond to with (without) antishadowing contributions. 
The dashed curves are the results of the DGLAP equation. 

Fig. 6 As Fig. 5 but for the sea-quark distribution function in proton. 

Fig. 7 The gluon recombination scale Qjiix) in the MD-DGLAP equation for proton 
(solid curve), Ca{A — 40) (dotted curve) and Pb(A = 208) (dashed curve); where the 
input distributions are Eqs. (3) and (4). 

Fig. 8 The predictions for the gluon distribution function in Ca(A = 40). The solid 
(dotted) curves are the results of the MD-DGLAP equation with (without) antishadowing 
contributions; The dashed curves are the results of the DGLAP equation; where the input 
distributions are Eqs. (3) and (4). 

Fig. 9 As Fig. 8 but for the sea-quark distribution function. 

Fig. 10 The dependence of the gluon distribution function in Ca(A = 40) at fixed 



values of a;, from with (solid curves) and without (dotted curves) antishadowing correc- 
tions in the present work. The dashed curves are the results of the DGLAP equation. 
Notice the plateau at small x] where the input distributions are Eqs. (3) and (4). . 

Fig. 11 As Fig. 10 but for the sea-quark distribution function. 

Fig. 12 As Fig. 8 but for Pb(^ = 208). 

Fig. 13 As Fig. 9 but for Pb(A = 208). 

Fig. 14 As Fig. 10 but for Pb(^ = 208). 

Fig. 15 As Fig. 11 but for Pb(A = 208). 

Fig. 16 Parts of Figs. 14 and 15. The hues show that the altitude of the plateau in 
C&{A = 40) and Ph{A = 208) almost grows linearly with increasing Q'^. 

Fig. 17 The predictions for the gluon distribution function in Ca{A — 40) and Ph{A — 
208) using the MD-GDLAP equation and input distributions Eqs. (10)-(13). 

Fig. 18 The predictions for the quark distribution functions in Ca,{A — 40) and 
Ph{A = 208) using the MD-GDLAP equation and input distributions Eqs. (10)-(13). 

Fig. 19 As Fig. 7 but the input distributions are taken as Eqs. (10)- (13). 

Fig. 20 The predictions for the gluon distribution function using the GLR-MQ evo- 
lution equation with a GRV-like input parton distributions. The solid and dotted curves 
are the results in proton and Pb(A = 208). 

Fig. 21 As Fig. 20 but for the sea-quark distribution function. 
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